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?3 i 1. Introduction predicts that, for a spherically symmetrical collapse, the 

radiation pressure produced by the massive YSO enor- 
A widely accepted, observationally supported, "standard" mous l um i nos ity would be sufficient to inhibit any further 
theory exists for the formation of low-mass, solar-type accret i n. That leads to the paradox that massive stars 
stars, involving inside-out collapse of molecular cloud frag- ( M > g Mq ) might never form, in disagreement with 
ments and formation of accretion disks (e.g. |Shu et alj the observed upper limits of hundreds of solar masses. 
U2S2I). For high-mass (> 10 M ) stars, on the contrary, In order to circumven t this dilemma, it has been pro- 
further effort (from both a theoretical and observational poged that high-mass stars might build up via merging 
point of view) is required to improve our understand- of lo er mass (M < 8 M fl ) YSOs (coalescence model 
ing of their formation process. A fundamental question _ e g iBonnell fc Batdr2002h . or, alternatively, by accret- 
is whether the mode of high-mass formation is a scaled- ing mater i a l from the infalling envelope via a thin disk, 
up version of the low-mass forming mode or weather it as j t occurs in low-mass stars (accretion model - e.g. 
involves different physical processes. This question arises lYorke fc Souunalterl2002h . In the disk scenario, in order to 
because stars with masses larger than 8 M Q contract to allow the gas of the disk to accre te onto the stellar surface, 
the main sequence and ignite while still accreting matter the removal of the excess angu l ar momentum is thought 
from the infalling envelope. Then, the "standard" theory to be " some how" associated to the ejection of a power- 

ful bipolar outflow, collimated along the disk-axis. Then 
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Abstract. We have conducted phase-reference multi-epoch observations of the 22.2 GHz water masers using the 
Very Long Baseline Array (VLBA) and multi-frequency study of the continuum emission using the Very Large 
Array (VLA) towards the high-mass star forming region (SFR) AFGL 5142. 29 maser features were identified and 
most of them were persistent over the four observing epochs, allowing absolute proper motions to be determined. 
The water maser emission comes from two elongated structures (indicated as Group I and Group II), with the 
measured proper motions aligned along the structures' elongation axes. Each group consists of two (blue- and 
red-shifted) clusters of features separated by a few hundreds and thousands of AU respectively for Group I and 
Group II. 

The maser features of Group II have both positions and velocities aligned along a direction close to the axis of 
the outflow traced by HCO + and SiO emission on angular scales of tens of arcsec. We predict that the maser 
emission arises from dense, shocked molecular clumps displaced along the axis of the molecular outflow. The two 
maser clusters of Group I are oriented on the sky along a direction forming a large angle (2; 60°) with the axis 
of the jet/outflow traced by Group II maser features. We have detected a compact (8.4 and 22 GHz) continuum 
source (previously reported at 4.9 and 8.4 GHz) that falls close to the centroid of Group I masers, indicating 
that the source ionizing the gas is also responsible for the excitation of the water masers. The kinematic analysis 
indicates that the Group I masers trace outflowing rather than rotating gas, discarding the Keplerian disk scenario 
proposed in a previous paper for Group I. Since the axis joining the two maser clusters of Group II does not cross 
the position of the continuum source, Group II masers might be excited by an (undetected) massive YSO, distinct 
from the one (pinpointed by the VLA continuum emission) responsible for the excitation of the Group I masers. 
Our results give support to models of accretion and jet ejection related to the formation of high-mass stars. 
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disk/outflow systems, whereas the alternative coalescence 
scenario would imply the destruction of ordered kinemat- 
ical structures surrounding the merging protostars. 

According to widely accepted theories, (proto)stellar 
jets may be launched by the combined action of mag- 
netic and centrifugal forces along the magnetic field 
lines threading the YSO's envelope. The magnetohydro- 
dynamical (MHD) models proposed to explain the ori- 
gin of (proto)stellar jets can be divided into two cate- 
gories, depending on the place where the jet is supposed 
to originate: 1) X-winds, which are magnetized stellar 
winds driven by the stell ar magnetic field , arising at the 
star/disk boundary (see IShu et ail l2000l for a review); 
2) disk-winds, where the jet stems from the surface of 
a Ke plerian disk threaded by the circumstellar magnetic 
field I Konigl fc Pudritzll2000l and references therein). 

A possible way, from an observational point of view, 
to distinguish between alternative (accretion and jet ejec- 
tion) models in massive star formation theories is study- 
ing the kinematics of the gas in the proximity (< 100 AU) 
of the YSO, corresponding to the "root" of the putative 
jet/disk system. An ideal diagnostic tool is provided by 
Very Long Baseline Interferometry (VLBI) observations 
at centimeter (cm) wavelengths of maser transitions of 
molecular species (such as OH, H 2 0, CH 3 OH). In par- 
ticular, for fast-moving 22.2 GHz water masers, multi- 
epoch VLBI observations, reaching angular resolutions of 
~0.5 mas (corresponding to ~1 AU at a distance of 2 kpc), 
permit accurate proper motions to be determined with 
time baselines as short as a few months. Combining the 
measured tangential velocities with the radial velocities 
derived via the Doppler effect, allows one to obtain the 
3D velocity distribution of the masing gas. 

Previous VLBI observations have shown that 22.2 GHz 
water masers are preferentially associated with colli- 
mated (jets) or wide-angle flows of gas at t he base of 
large r-scale molecul ar outflows (W3 IRS 5, llmai et all 
l2000t W75N-VLA1. iTorrelles etall l2003t S h 2-255 IR 
IRAS 23139+5939, WB 89-234. and OMC2. lOoddi et all 
120051 IRAS 20126+4104. lMoscadelli et al.ll2005l) . In a few 
objects CNGC 207 1 IRS 1 /3. ISeth et a,l.l2002t AFGL 5142, 
iGoddi et alJ2 004 - hereafter Paper I), accretion disks have 
been proposed to explain linear clusters (size ^10-100 AU) 
of H2O masers detected through VLBI observations. 

The 22.2 GHz H 2 masers in AFGL 5142 (for a de- 
tailed description of the source see Paper I) were observed 
with the European VLBI Network (EVN) at four epochs 
(from October 1996 to November 1997) and found to be 
distributed across two structures (indicated as "Group I" 
and "Group II"). Such structures have been interpreted 
as being possibly associated with a disk/jet system, in- 
dicating that AFGL 5142 may be an ideal candidate to 
test the formation theories of high-mass stars. However, 
the proposed interpretation needs more stringent observa- 
tional constraints, the multiepoch EVN observations suf- 
fering two major drawbacks: 1) too long time separation 
between consecutive epochs (> 3 months, more than the 
average life time of the maser features); 2) limited sen- 



sitivity (average detection threshold of ^0.3 Jy beam -1 ). 
These limitations prevented us to detect and measure the 
proper motions of the weakest and less-longeval maser fea- 
tures. 

Recently, we performed follow-up observations of the 
22 GHz H 2 masers in AFGL 5142 using the Very Long 
Baseline Array (VLBA), taking advantage of both high 
sensitivity and an optimized time separation between con- 
secutive epochs. The VLBA observations were performed 
in phase-reference mode to derive the absolute positions 
and proper motions of the water maser features. 

We have also performed multi-frequency VLA-B con- 
tinuum observations of AFGL 5142, aiming to precisely lo- 
cate the position of the exciting YSO(s), to put constraints 
on the nature of the previously observed radio contin- 
uum e mission llCarnepter e^ al.lll99fj iMcCutcheon^et alJ 

il. 



19911: ITorrelles et al l h 9921 iHi inter et alJ ll 9951 ll QQflt 

Carral et al. 1999), and, finally, to investigate the relation- 
ship between the water masers and the continuum source. 

Section describes our VLA and multi-epoch VLBA 
observations and gives technical details on the data analy- 
sis. Section [3] compares our results with previous interfer- 
ometrical observations of the region. In Section plausi- 
ble kinematic models for interpreting the measured posi- 
tions and velocities of the maser features are investigated. 
Conclusions are drawn in Section [S] 

2. Observations and data analysis 

2.1. VLA Observations 

The continuum observations were performed using the 
VLA in its B configuration on March 1 and 11 2005, for 
a total observing time of 4 hours. We observed a band- 
width of 50 MHz centered at a central frequency of 8.4601, 
22.4601, and 43.3399 GHz (V, K, and Q-bands, respec- 
tively). The data were edited, calibrated, and imaged us- 
ing the NRAO's Astronomical Image Processing System 
(AIPS) package. Absolute amplitude calibration was ob- 
tained from observations of 3C147 (with flux density of 
4.8, 1.8, and 0.9 Jy for X, K, and Q-band, respectively). 
In each band, the observations were performed using the 
"fast switching" method, where short (420, 180, 60 s) 
scans of AFGL 5142 were alternated with 60 s scans of 
the nearby phase calibrator J05181+33062. The error in 
the absolute calibration is estimated to be within 10% in 
X, 20% in K and 30% in Q-band. Table □ reports the 
main observational parameters (restoring beam FWHM 
sizes and position angle, and the map RMS noise level) 
for the three observing bands. 

2.2. VLBA Observations 
2.2.1. VLBA runs 

The observations in the 616 — 5 2 3 H 2 maser line at 
22.2 GHz were conducted using the VLBA at four epochs 
(16 October and 22 November 2003, 1 January and 8 
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Table 1. Radio continuum observations. 



Band 


Beam size 


Beam P. A. 


a 


q(J2000) 


(5(J2000) 


S v 


FWHM Size 


PA. 




(") 


(°) 


(/ijy beam -1 ) 


(h m s) 


(° ' ") 


(mJy) 


(arcsec) 


(°) 


X 


0.88 x 0.82 


71 


26 


05 30 48.019 


33 47 54.6 


0.75 ± 0.08° 






K 


0.38 x 0.38 





52 


05 30 48.020 


33 47 54.5 


0.7 ±0.1 


0.35 x 0.27 6 


130 


Q 


0.17 x 0.15 


31 


270 






< 0.9 C 







Note.- For each observing band, Cols. 2 and 3 give respectively the FWHM sizes along the major and minor axes, and the 
position angle of the restoring beam; Col. 4 the RMS noise level on the maps; Cols. 5 and 6 the R.A. and DEC absolute 
position of the intensity peak; Col. 7 the integrated flux density of the maps; Cols. 8 and 9 the deconvolved FWHM sizes and 
the position angle of the detected sources. 

a The absolute flux density uncertainty is estimated to be < 10% at A-band, < 20% at if-band, < 30% at Q-band. 

6 The source is compact at all observed frequencies, with the exception of the Tf-band, for which the deconvolving routine 

(SAD) produced FWHM sizes comparable with the observing beam. 

c At 43.3 GHz the flux density upper limit corresponds to 3a. 



February 2004), each epoch lasting for 12 hours. The ob- 
servations were performed in phase-reference mode, al- 
ternating scans on the maser source, AFGL 5142, and 
the phase-reference source, J0518+3306, with a switch- 
ing cycle of 70 s. J05181+33062 is separated 2.6° from 
AFGL 5142, belongs to the VLBA calibrator catalog, 
and has very accurate coordinates (R.A. and DEC un- 
certainties < 1 mas). Interlapsed every «80 min, 3-min 
scans on several continuum sources (0528+134, 0552+398, 
0333+321, 0642+449) were observed for calibration pur- 
poses. The total on-source integration time for the maser 
target was about 5.5 hours. 

Both circular polarizations were recorded using a 
16 MHz bandwidth centered on the LSR velocity of - 
4.8 km s _1 . The data were correlated with the VLBA FX 
correlator in Socorro (New Mexico) with an integration 
period of 1 s. The correlator used 1024 spectral channels 
corresponding to a channel separation of 0.2 km s _1 . 

2.2.2. Data reduction 

Data reduction was performed using the AIPS package, 
following the standard procedure for VLBI line data. Total 
power spectra of the continuum calibrators were used to 
derive the bandpass response of each antenna. The am- 
plitude calibration was performed using the information 
on the system temperature and the gain curve of each 
antenna. 

For each observing epoch, a single scan of a strong 
calibrator was used to derive the instrumental (time- 
independent) single-band delay and the phase offset be- 
tween the two polarizations. After removing the instru- 
mental errors, all calibrator scans were fringe-fitted to de- 
termine the residual (time-dependent) delay and the fringe 
rate. The corrections derived from calibrators were applied 
to the strongest maser component, chosen to refer the vis- 
ibilities of all the other maser emission channels. In each 
observing epoch the strongest maser component emitted 
at the same velocity, Vlsr = -4.9 km s _1 , and exhibited a 
simple spatial structure consisting of a single, almost un- 
resolved spot. The visibilities of the reference channel were 



fringe-fitted to find the residual fringe rate produced both 
by differences in atmospheric fluctuations between the cal- 
ibrators and the maser, and by errors in the model used at 
the correlator. After correcting for the residual fringe rate, 
the visibilities of the reference channel were self-calibrated 
to remove any possible effect of spatial structure. Finally, 
the corrections derived from the reference channel were 
applied to data of all spectral channels. 

Spectral channel maps were produced extending over a 
sky area of (Aa cosS x AS) 6" x 6" and covering the whole 
velocity range where signal was visible in the total-power 
spectra (from —10 to 5 km s -1 ). At each epoch, the maser 
emission centres are found to be distributed within an area 
of 2"xl."5. The CLEAN beam was an elliptical Gaussian 
with a FWHM size of 0.7 x 0.4 mas. In each observing 
epoch, the RMS noise level on the channel maps, a, varied 
over a similar range of values, 4-40 mJy beam -1 . 

2.2.3. Absolute position determination 

Auto-referencing the visibilities phases of the maser data 
causes the loss of the absolute position information. Since 
in our case the water maser source is significantly stronger 
than the phase reference source, J0518+3306, we first 
attempted to determine the maser absolute position us- 
ing the reverse phase referencing technique. All the phase 
calibration (fringe-fitting + self-calibration) solutions ob- 
tained working on the reference maser component were 
applied to the J05181+33062 data. So doing, the position 
of the phase-reference source is shifted by a vector equal in 
absolute value but with opposite sign respect to the offset 
of the maser reference spot from its correlated position. 
Finally, the maser-referenced visibilities of J05181+33062 
were Fourier-transformed, but the resulting map produced 
only a marginal detection (with a peak flux density of 
- 10 mJy). 

Then, we experimented the alternative way of 
direct phase-referencing, fringe-fitting the scans of 
J05181+33062 and applying the obtained phase calibra- 
tion solutions to the maser data. The image constructed by 
Fourier transforming the visibilities of the maser source, 
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Table 2. Offsets of the reference feature from the correlated position: a(J2000) = 05 h 30 m 48."033, <5(J2000) = 
33°47'54."6; Col. 1 reports the observing epoch; Cols. 2 and 3 respectively the R.A. and DEC apparent offsets, as 
derived from the phase-reference calibrated maps; Cols. 4 and 5 the corrections for the combined action of Galactic, 
Sun and earth motion; Cols. 6 and 7 the corrected (R.A. and DEC) offsets. 



Epoch 


Aai 


A5i 


Aa 2 


A5 2 


Aatf 


A5 f 




(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


(mas) 


16 Oct 03 


-83.1 (0.4) 


-244.7 (0.2) 


-1.2 


1.0 


-84.3 (0.4) 


-243.7 (0.2) 


22 Nov 03 


-83.4 (0.1) 


-244.6 (0.1) 


-0.9 


1.2 


-84.3 (0.1) 


-243.4 (0.1) 


1 Jan 04 


-83.7 (0.3) 


-245.3 (0.2) 


-0.6 


1.4 


-84.3 (0.3) 


-243.9 (0.2) 


8 Feb 04 


-84.0 (0.3) 


-245.5 (0.2) 


-0.3 


1.7 


-84.3 (0.3) 


-243.8 (0.2) 



calibrated with the phase-reference continuum solutions, 
presented a signal-to-noise ratio sufficiently high to allow 
us to determine the offset of the reference maser (from 
the position used at the correlator) with an accuracy of 
tenths of milliarcsec. For the first observing epoch, the 
derived positional offset is: Aai = -0."0831 ± 0."00035, 
A6 1 = -0."24469 ± 0."0002. The uncertainty on the posi- 
tional offset is estimated by the ratio of the FWHM beam 
size to the dynamic range of the map. In order to obtain 
the peculiar motion of the maser source, the offsets derived 
at each epoch have been corrected for systematic effects 
associated to Galactic rotation, annual parallax and solar 
motion. The corrected offsets of the reference feature are 
reported in Table El 

2.2.4. Identification of maser emission 

Each channel map of the maser source was searched for 
emission above a detection threshold taken equal to the 
absolute value of the minimum in the map (varying across 
the channels in the range 5-7 a) . The detected maser spots 
were fitted with two-dimensional elliptical Gaussians, de- 
termining position, flux density, and FWHM size of the 
emission. Consistently to Paper I, the term of "feature" is 
used to indicate a collection of spectrally and spatially 
contiguous maser spots. A maser feature is considered 
real if it is detected in at least two contiguous channels 
(i.e., line-width > 0.4 km s -1 ), with a position shift of 
the intensity peak from channel to channel smaller than 
the FWHM size. The VLBA observations provided an im- 
provement in sensitivity of about an order of magnitude 
respect to the previous EVN observations, lowering the 
detection threshold to ~ 30 mJy beam -1 and allowing 
the detection of a much larger number of maser spots 
and features. In each observing epoch, about 270 maser 
spots were identified, grouped into 29 distinct maser fea- 
tures. Table |3 gives the parameters of the identified fea- 
tures. The absolute (R.A. and DEC) positional offsets of 
a given feature, evaluated with respect to the maser po- 
sition used at the correlator (a(J2000) = 05 /l 30 m 48."033, 
<5(J2000) = 33°47'54."6), are calculated with a two step 
procedure. First, relative positions (respect to the refer- 
ence maser spot) are derived from the (error-weighted) 
mean of the maser spot positions contributing to the fea- 
ture emission. Then, at each epoch the relative positions 



are combined with the absolute positional offsets of the 
reference spot (reported in Table[5J). For each feature, the 
absolute positional errors are estimated taking the sum 
of the feature relative positional uncertainties (evaluated 
by the weighted standard deviation of the spot positions) 
with the absolute positional uncertainties of the reference 
spot. 

2.2.5. Proper motion measurements 

For the features persistent over three or four epochs, the 
proper motions have been calculated performing a (error- 
weighted) linear least-squares fit of the positional offsets 
with time. For the VLBA observations the time separation 
between consecutive epochs was shorter than for the pre- 
vious EVN observations (~1 month vs 3-6 months) and 
that allowed to measure the proper motions of the most 
variable features, increasing significantly the number of 
measured proper motions (23 vs 2 of previous EVN ob- 
servations, considering only features persistent over > 3 
epochs). For a representative sample of persistent fea- 
tures, Fig. ^ shows the time variation of the R.A. and 
DEC (relative) offsets and the best linear fit giving the 
proper motion. Among the features observed at three or 
more epochs, the proper motions are derived only for those 
moving in a straight line at constant velocity (within the 
positional errors). Adopting such a criterion, no proper 
motion is derived for the feature with label number " 26" , 
although detected at all the four observing epochs. 

The derived absolute proper motions are reported in 
Table 01 The uncertainties of the absolute proper motions 
are the formal errors of the linear least-squares fits and 
are generally much larger than those associated with the 
relative proper motions, owing to the relatively large error 
of the absolute position of the reference spot. The phase- 
reference calibrator is too weak and too much detached 
from the maser target to allow an effective removal of the 
atmospheric noise affecting the maser visibilities. 

3. Observational results 

3.1. Radio continuum emission 

The VLA-B, X and if-band continuum images of 
AFGL 5142, above a detection threshold of 5er, show 
a single, compact source coinciding (within the posi- 
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Fig. 1. Measured relative proper motions of the VLBA H 2 maser features in AFGL 5142. For each of the time 
persistent features, the top and the bottom panels report the time variation of, respectively, the DEC and the R.A. 
offsets (relative to the reference feature "1"). In each panel, the dotted line shows the proper motion calculated by the 
(error-weighted) linear least-squares fit of the positional offsets with time. No proper motion is derived for the feature 
with label number "26", since it does not show a linear variation of positional offsets with time. 
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Table 3. Parameters of maser features detected with the VLBA in AFGL 5142. 





Feature 


Vlsr 


Fint 


Aa 


AS 


Vi 


Vy 


Vmod 






(km s _1 ) 


(Jy) 


(mas) 


(mas) 


(km s _1 ) 


(km s" 1 ) 


((km s" 1 ) 




1 


-5.0 


33.2 


-84.3 (0.4) 


-243.7 (0.2) 


0.8 (11) 


-11 (7) 


11 (7) 




2 


-5.7 


7.1 


-59.6 (0.4) 


-242.9 (0.2) 


14 (14) 


-15 (7) 


20 (11) 




3 


-6.3 


5.3 


-191.1 (0.4) 


-279.2 (0.2) 


-9 (13) 


16 (7) 


18 (9) 




4 


-7.8 


1.1 


-66.6 (0.4) 


-237.8 (0.2) 


5(14) 


-11 (7) 


12 (9) 




5 


-4.4 


0.1 


-57.9 (0.3) 


-242.1 (0.2) 










6 


-3.1 


0.6 


-248.1 (0.4) 


105.5 (0.2) 


-9 (13) 


8(7) 


12 (11) 




7 


-2.7 


1.1 


-77.8 (0.4) 


242.3 (0.2) 


8(13) 


-10 (7) 


13 (10) 




8 


0.1 


0.5 


-234.7 (0.4) 


0.0 (0.2) 


-9 (13) 


2 (7) 


9 (13) 




9 


-0.2 


0.6 


-233.2 (0.4) 


-0.4 (0.2) 


-6 (13) 


4(7) 


8 (12) 


Group I 


12 


-1.4 


0.5 


-186.7 (0.4) 


-26.2 (0.2) 


-26 (16) 


9(8) 


28 (16) 




13 


-2.2 


0.3 


-182.8 (0.4) 


-27.4 (0.2) 


0.9 (13) 


9(7) 


9(7) 




14 


-1.3 


5.2 


-299.9 (0.4) 


43.5 (0.2) 


-13 (13) 


5 (7) 


14 (13) 




15 


-1.7 


0.7 


-253.9 (0.4) 


92.1 (0.2) 


-7 (13) 


6(7) 


9(11) 




16 


-1.9 


0.3 


-290.0 (0.4) 


-24.0 (0.2) 










17 


-1.7 


0.5 


-281.9 (0.4) 


-99.0 (0.2) 


-6 (14) 


-5 (7) 


8(11) 




18 


4.6 


0.9 


-171.4 (0.4) 


-169.9 (0.2) 


-18 (14) 


4(7) 


18 (13) 




i n 


1 . 9 


n r 
U.O 


—222.4 (O.OJ 


4.U (U.Z) 


-0.3 (16) 


11 (8) 


11 (8) 




20 


-2.0 


0.1 


-303.2 (0.4) 


7.0 (0.5) 


-12 (13) 


-2 (20) 


13 (14) 




21 


-1.3 


0.2 


-184.5 (0.4) 


-26 (0.2) 


-18 (13) 


16 (7) 


24 (11) 




22 





0.3 


-291.0 (0.2) 


-5.8 (0.1) 


-26 (17) 


5 (9) 


27 (17) 




23 


1.8 


0.1 


-290.1 (0.3) 


-19.96 (0.2) 










24 


1.3 


0.1 


-241.4 (0.2) 


118.3 (0.1) 


-5.0 (17) 


3(9) 


6 (15) 




25 


-0.6 


0.1 


-240.6 (0.3) 


118.9 (0.2) 










10 


-0.1 


0.7 


199.8 (0.4) 


89.2 (0.2) 


22 (14) 


11 (7) 


24 (13) 




11 


0.9 


0.1 


197.4 (0.4) 


90.6 (0.2) 


6 (13) 


19 (7) 


20 (8) 


Group II 


26 


-7.4 


6.6 


-469.2 (0.4) 


-1587.9 (0.2) 










27 


4.1 


0.5 


-450.0 (0.2) 


-1481.0 (0.1) 










28 


-5.7 


4.9 


-921.5 (0.4) 


-1813.4 (0.2) 


-17 (14) 


-16 (7) 


24 (11) 




29 


-4.2 


0.7 


-921.6 (0.4) 


-1812.9 (0.2) 


-27 (13) 


-18 (7) 


32 (12) 



Note.- For each identified feature, Col. 1 gives the label number; Cols. 2 and 3 the line-of-sight velocity and the integrated 
flux density of the highest-intensity channel; Cols. 4 and 5 the absolute positional (R.A. and DEC) offsets (and, within 
brackets, the associated errors), evaluated with respect to the maser position used at the correlator (a(J2000) = 05 ,I 30 m 48. 
"033, <5(J2000) = 33°47'54."6); Cols. 6, 7 and 8 the projected components along the RA and DEC axes and the absolute 
value of the derived proper motions (and, within brackets, the associated errors). 



tional er rors) with the previously detected emission at 



4.9 GHz dCaroenter et alJll990: McCutcheon et all 


1991 


and 8.4 GHz 


Torrclles et al. 1992; Hunter et al. 


1995 


1999; Carral et al. 


1999). No emission was detected at 



43 GHz and only a 3cr upper limit to the source flux den- 
sity could be derived. Table gives the measured posi- 
tions and flux densities for the three observing bands. At 
K-band, the emission appears slightly resolved, with a de- 
convolved elliptical Gaussian shape having a FWHM size 
(estimated using the AIPS task SAD) of 0."35x0."27 (cor- 
responding to a diameter of ~ 600 AU), and elongated 
along a southeast-northwest (SE-NW) direction. 

Previous observations at 8.4 GHz continuum, con- 
ducted using the VLA in A llHunter et all 1*1 995). B 
1 Hunter et all 1 19991) . C (|Torrelles et alJ Il992|) . and D 
I Carral et a 1. 1 119991) configurations, all s how a compact 
continuum emission, with the exception of lTorrelles et alJ 



(1992), who presented marginal indication for the emis- 
sion being extended at scales of ~ 20" 1 (see Fig. 3 of 
iTorrelles et al 1ll992l) . Fig.Hreports our and previous con- 
tinuum observations towards AFGL 5142. At 8.4 GHz we 
have reported only our VLA -B measurement, not consid- 
ering the previous value of iHunter et al.l l)l999l) . which 
is much more uncertain 2 . One notes that, at both 4.9 
and 8.4 GHz, the integrated flux density increases with 
the FWHM beam size of the observation, suggesting the 
presence of an extended component with size > 4" (the 

1 In the lTorrelles et alJ il992T) map, the presence of extended 
emission is suggested only by the weakest contours, below the 
3<r level. 

2 Th e 8.4 GHz observations were conducted bv lHunter et alJ 
( 1999J) in snapshot mode, with a few minutes of integration, 
whereas our 8.4 GHz observations last for about one hour, 
allowing a significant sensitivity improvement. 
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Fig. 2. The radio spectral energy distribution of 
AFGL 5142. Filled triangles are our VLA-B measure- 
ments of the flux density at 8.4 (0.75 mJy with a rj 0."9 
beam), 22.4 (0.7 mJy with a « 0."4 beam), and 43.3 GHz 
(< 0.9 mJy with a ss 0."2 beam; the upper limit is in- 
dicated by the downward arrow). The asterisks are pre- 
vious flux density measurements at 4.9 GHz (1.0 mJy 
with a m 12."5 beam using VLA-D - ICarpenter et alfl 990 
(C90) - and 0.5 mJy wi th a w 4" beam using VLA-C - 
iMcCutcheon et alJll99ll (M91)) and 8.4 GHz (0.83 mJy 
with a ss 0."3 beam using VLA-A - iHunter etaD Il995l 
(H95); 1.0 mJy wit h a « 7."7 beam using VLA-C - 
iTorrelles et"al]ll992l (T92): 1.2 m .Tv with a w 7"5 beam 
using VLA-D - ICarral et aflll999l (C99)). The size of the 
symbols scales logarithmically with the FWHM size of the 
beam. 



VLA-C beam). Going from 4.9 to 8.4 GHz the VLA-C 
flux density increases by a factor of 2 (less than a factor 
of 3, as expected for optically thick emission), whereas 
the VLA-D flux density remains substantially constant. 
That is consistent with the presence of a compact core 
of emission at 4.9 GHz, that becomes optically thin at 
8.4 GHz. Accordingly, our VLA-B measurements indi- 
cate that the integrated flux density remains constant 
or slightly decreases at frequencies larger than 8.4 GHz. 
Assuming an average optical depth w 1 at 8.4 GHz for the 
emission inside the VLA-B beam («0."9), one finds that 
thebeam- averaged value for the density of the ionized gas 
is 2 x 10 5 cm~ 3 . 

Basing on the extended nature of the radio continuum 
and the derived high value for the density of the ionized 
gas, the most likely interpretation is that of a young HII 
region. To account for the measured radio flux density, 
an ionising flux by a ZAMS star with spectral type B2 
or earlier is required, corre sponding to an expected value 
of ste llar mass > 10 M© l(Vacca et al.lll99fil Ipalla et alJ 



3.2. 22.2 GHz water maser emission 

Fig-EDcompares our VLA and VLBA results with previous 
interferometric observations. Top right-hand panels show 
the high-velocity molecular outflows seen in HCO + (1 — > 
0) an d SiO (v = 0, 2 -> 1) with OVRO iHunter et alJ 
Il999h . The 22.2 GHz water masers in AFGL 5142 (prior 
of our EVN obser vations) were observ ed wit h the VLA at 
two e pochs (!992. lHunter et alJll995 1998. lHunter et alJ 
119991) . and found within a few arcseconds from a compact, 
thermal continuum source (observed both with the VLA 
at 8.4 GHz and with OVRO at 88 GHz), located at the 
center of the SiO-HCO + outflow (Fig. upper left-hand 
panel). The lower panel of Fig. 01 shows the positions and 
the velocities of the H2O maser features as derived by 
our multi-epoch VLBA observations, overlaid on our VLA 
map of the 22 GHz continuum emission. 

Basing on the spatial distribution and the proper mo- 
tion orientation of the detected maser features, the water 
maser emission appears to trace two elongated structures, 
oriented at quite different P.A. (AP.A.> 60°). Using the 
same nomenclature of Paper I, we refer to such structures 
as Group I and Group II. In Table we report the de- 
tected maser features separately for each Group. 

Masers of Group I are found closer (within «300 mas) 
to the 22 GHz continuum emission peak. It is worth noting 
that, basing on the accurate absolute positions we have 
determined for the water masers and the thermal radio 
source, the 22 GHz continuum peak falls right at the center 
of the maser distribution of Group I (Fig. |21 lower panel) . 
The features detected towards SE of the continuum source 
have LSR velocities blue-shifted (respect to the LSR ve- 
locity of the region) whereas those detected towards NW 
are red-shifted. The measured transverse velocities of the 
blue and red-shifted features have similar amplitudes (in 
the range 6-27 km s _1 ) but opposite versus of motion. 

Maser features of Group II are concentrated in two 
clusters, having a sky-projected distance of ~ 1700 mas 
(corresponding to ~ 3500 AU at a distance of 1.8 kpc). 
The line connecting the two clusters is oriented along 
a northeast-southwest (NE-SW) direction and the aver- 
age line-of-sight velocity of the features detected towards 
northeast is red-shifted whereas that of the southwestern 
features is blue-shifted. The measured proper motions in- 
dicate that the two clusters are moving away from each 
other along a NE-SW direction with similar velocities, in 
the range 20-32 km s _1 . 

Looking at Fig. (bottom panel), one can note that 
the line joining the two maser clusters of Group II does 
not cross the position of VLA 22 GHz continuum peak, 
but it is instead shifted towards S-E by >0."4. 



4. Discussion 

The comparison of the spatial and velocity distribution of 
the molecular gas traced by the HCO + and SiO emission 
and the Group II water masers, suggests that these latter 
may be tracing the inner portion of the large-scale molec- 
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Fig. 3. VLA and VLBA results towar ds AFGL 5142. ( Up per right-h and panel) Contour maps of the OVRO high-velocity emission of HCO+ (1 -» 
0) (left) and SiO (v = 0, 2 — > 1) (right) IHunter et al.H999fl . ( Upper left-hand panel) Positions of the 22 GHz masers dete cted with the VLA i n 199 2 
(open squares) and 1998 (filled dots); OVRO 88 GHz (boxed cross) and VLA 8.4 GHz (plain cross) continuum sources IHunter et aH199afT999l) . 
(Bottom panel) Positions of the VLBA H2O mascr features are marked with filled triangles, different colours denote the features' LSR velocities, 
according to the colour scale on the right-hand side of the panel, and the arrows indicate the measured absolute proper motions, whose amplitude 
scale is given at the bottom of the panel; the labels lb, Ir and lib, Hr individuate the blue-shifted and red-shifted mascr clusters of Group I and 
Group II, respectively; the contour map (with contours representing 5, 6, 7, 8, and 9 times the 50 fiJy beam -1 RMS noise) shows the 22 GHz VLA 
continuum emission; the filled circle in the upper right corner represents the restoring beam of the VLA-B 22 GHz continuum map. 
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ular outflow. This interpretation finds support upon three 
arguments: the elongated distribution of the masers is al- 
most parallel to the position angle of the HCO + outflow; 
the blue-shifted features lie northward and the red-shifted 
ones lie southward, qualitatively matching the velocity dis- 
tribution of the large-scale molecular outflow; the orien- 
tation of the measured proper motions is approximately 
aligned with the axis of the molecular outflow. 

In Paper I, basing on low-sensitive EVN data, the ob- 
served velocity pattern of Group I masers has been inter- 
preted in terms of a Keplerian rotation. In order to test the 
Keplerian hypothesis, we fitted the positions and the 3D 
absolute velocities of maser features with the Keplerian 
disk model described in Paper I (see Appendix ^] for a 
more detailed description of the model formalism). The 
minimizing algorithm (Eq. IA.1I in App. [SJ converged to 
a best fit solution. The geometrical parameters of the 
best fit disk are found to be in good agreement with 
those previously derived with the EVN maser data (see 
Paper I). Since the measured relative velocities are in our 
case three-five times more accurate than absolute veloci- 
ties (see Sect. l2.23fr . we repeated the fit comparing model 
predicted and measured relative (to the reference maser 
feature) velocities fEa. IA.21 in App. 0. In this case, the 
Keplerian model failed in finding a minimum, indicating 
that the measured velocity field for the Group I masers 
is actually incompatible with Keplerian rotation. This re- 
sult demonstrates that the solution obtained fitting both 
absolute (VLBA) and relative (EVN) proper motions is a 
spurious result owing to the large measurement errors. 

The unsuccess of the Keplerian model is not surpris- 
ing. Looking at the Fig. one notes that the motion 
traced by Group I water masers appears to be expansion 
from (rather than rotation around) the putative location 
of the YSO, as marked by the peak of the thermal con- 
tinuum emission. In particular, proper motions in both 
Group I (SE and NW) maser clusters are concordly ori- 
ented outward from the expected YSO position, whereas 
in a Keplerian rotation they should be directed both in- 
ward and outward, unpreferably tracing both the far and 
near side of a disk annulus. 

To test if Group I masers may trace an outflowing mo- 
tion, the data were fit ted with the conica l outflow model 
described in detail in iGoddi et al.l (|2005l) . Since the ab- 
solute values of the measured velocities do not present a 
large variation (compared with the observational errors - 
see Table 0), in the conical outflow model we assumed the 
maser features to move with a constant vel ocity (differ- 
ently from the Hubble motion assumption in lGoddi et alJ 
120051) . The best fit solution was looked for minimizing the 
\ 2 given by the sum of the square differences between the 
model and measured, relative, 3D velocities (Eq. I A. 21 in 
App. 0). The fit parameters were searched over a wide 
range of physically plausible values: the cone vertex (i.e., 
the YSO location) within a distance of 300 AU from the 
22 GHz continuum peak; the cone axis orientation over 
the full 4-7T solid angle; the cone aperture within the in- 
terval 10 — 85°; the constant expansion velocity from 1 



to 50 km s _1 . The fit parameters were varied with steps 
comparable or smaller than the associated fit errors. The 
minimization algorithm converged towards the best fit so- 
lution, whose main characteristics are hereafter briefly re- 
sumed. The best fit conical jet is found to have a large 
(« 80°) opening angle. The best fit cone-axis is predicted 
to form an angle of ~ 30° with the line-of-sight and to 
have a position angle projected onto the plane of the sky 
of w 28°. The (constant) expansion velocity has an abso- 
lute value, vq, of 13 km s _1 . 

At this point a fundamental question arises: what 
might be the driving source of the observed flow motions 
of Group I and II masers? 

First we consider the scenario where an unique YSO 
is responsible for the radio continuum emission, the ex- 
citation of both maser groups, and the acceleration of 
the large scale molecular outflow. Si nce Group I individ- 
uates a dense (riHo > 10 7 cm -3 - lElitzur et al.1 Il989t 
IKaufman fc Neufeldlll996h clump of gas around the ex- 
pected YSO position, oriented on the sky at large an- 
gle from the axis of the jet/outflow system traced by 
the Group II maser and (at larger scale) by the HCO + - 
SiO emission, it is presented the possibility that Group 
I masers may originate in the accretion disk surround- 
ing the YSO. We speculate that Group I maser features 
might emerge from the atmosphere of a flared disk (with 
a disk scale height H(r) scaling linearly with r), in the 
region where modern MHD models predict the launch of 
the protostellar jet. As stressed in Sect. 1, such models can 
be divided into two categories, X-wind and disk-wind. In 
the former model, the jet is expected to originate within 
0.1 AU from the YSO and to collimate at radii < 10 AU, 
whereas the latter model predicts the wind to be launched 
across a significant range of disk radii (from tens of AU 
up to several hundreds of AU, for high-mass YSOs) and 
to collimate along the disk axis at correspondingly larger 
distances from the YSO. The detected H2O masers are dis- 
tributed at distances from the YSO position in the range 
10-600 AU. These large distances are clearly incompat- 
ible with the X-wind requirements, but they might still 
find an explanation within the disk-wind scenario. In this 
view, the disk- wind model would explain both Group I and 
Group II water masers, these latter tracing the collimated 
portion of the jet-outflow system at larger distances from 
the YSO (> 1000 AU). 

In the proposed scenario of a single powering source, 
one would expect the NE-SW outflow axis to be centered 
in between the maser clusters Ir and lb, i.e. at the position 
of the continuum source. That is not the case and this is 
a clear difficulty for the disk-wind interpretation. 

The alternative scenario is the one where the water 
masers in AFGL 5142 trace two physically uncorrelated 
jets (one directed towards SW-NE and the other one 
towards SE-NW), each jet being powered by a distinct 
YSO. As already discussed in Sect.EU the YSO respon- 
sible for the excitation of both the ionized source and the 
Group I masers is expected to be massive (> 10 M@). 
The other YSO should power the Group II masers and 
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the HCO + -SiO outfl ow, which contains a mass of 35 M Q 
ijHunter et aljfl99!^ . From excitation and energetic con- 
siderations one would expect this YSO to be also a mas- 
sive one. Geometrical considerations also suggest that this 
YSO may be located on the axis of the outflow traced by 
the Group II masers. Along this axis, separated by more 
than 0."4 of the 22 GHz continuum peak, no compact con- 
tinuum emission was detected by our sensitive VLA obser- 
vations. A possible explanation is that the YSO exciting 
Group II masers might be found in an evolutionary stage 
earlier than the YSO powering the Group I masers, when 
the ionized gas is still too much confined by the dense 
circumstellar environment to be detectable. 

We have reviewed the literature looking for indica- 
tions of a multiple system of YSOs in the r egion of de- 
tection of th e water masers. IHunter et all l|l995|) and, 
very recently. Ichen et aD l|2005|) performed NIR JHK and 
H2 V = 1 — 5(1) imaging observations of the molecu- 
lar cloud core in AFGL 5142, revealing the presence of a 
cluster of YSOs distributed over a region of diameter of a 
few arcminutes. H2O masers ar e found in correspo ndence 
of the infrared so urce NIRS1 oflChen et alJ l|2005|) (IRS1 
in the notation of lrlunter et al.lll995|) . The minimum dis- 
tance between NIRS1 a nd another memb er of the cluster 
is > 5" (see Table 2 of IChen et aDl2005l) . whereas (con- 
sidering the diameter of the maser distribution) the two 
(putative) YSOs powering the two maser groups should 
be separated by a distance 0."4 ;$ d 2" (corresponding 
to 7005; d < 3600 AU). Hence, the available NIR observa- 
tions, having an angular resolution > 1", do not allow us 
to identify two distinct powering sources for maser Groups 
I and II. Towards AFGL 5142, considering an angular 
scale larger than that presented in Fig. [21 (top right-hand 
panels), a complex pattern of molecular outflows has also 
been found, which might sugges t the p resence of a mul- 
tiple YSO system: IHunter et aD l| 19951) detected two CO 
outflows, nearly perpendicular to each other, one, directed 
towards SE-NW, extended on scales of a few arcminutes, 
the other one, direc ted towards SW-NE, more compact; 
IHunter et aD <|l999j) detected a HCO+-SiO outflow elon- 
gated in a N-S dir ection (Fig. 03 upper right-hand panels); 
IChen et al.l §005) found H2 jet-li ke structure s along each 
of the three molecular outflows. IChen et al.l l|2005|) pro- 
pose that the SW-NE oriented jet/outflow system is pow- 
ered by NIRS1, whereas the SE-NW oriented jet/outflow 
system is driven by NIRS30, 30" apart from NIRS1 (cor- 
responding to the IRAS source in the AFGL 5142 SFR). 
Since NIRS1 and NIRS30 are by now the two u nique recog- 
nized massive YSOs of the AFGL 5142 cluster l|Chen et alJ 
2005), the available molecular observations (with an an- 
gular resolution > 1") provide no hints for the presence of 
two massive YSOs within the region of a few arcsec where 
the water maser emission is observed. 

In conclusion, even if the geometry of the outflows 
traced by the two detected maser Groups may favour 
the interpretation of two distinct exciting YSOs, the data 
presently available do not permit unambiguously to iden- 
tify more than a single YSO in the region of maser emis- 



sion. In this respect, the axes of the outflows traced by 
water masers on scales of hundreds and thousands of AU 
can be good places where to look for the YSOs and the 
present water maser observations can drive future studies 
of the SFR AFGL 5142. 

5. Conclusions 

This article reports the results of VLBA multi-epoch ob- 
servations of the 22.2 GHz water masers and VLA multi- 
frequency study of the continuum emission towards the 
SFR AFGL 5142. Respect to the previous EVN obser- 
vations, the VLBA provided a significant improvement 
in sensitivity (by a factor ~ 10) and an optimized time 
separation between consecutive epochs (~1 month vs 3-6 
months), resulting in a larger number of detected maser 
features (29 vs 12) and measured proper motions (over 
three or four epochs, 23 vs 2). The VLBA observations 
were performed in phase-reference mode, allowing the de- 
termination of absolute positions and proper motions of 
the water maser features. 

The water maser emission is found to originate from 
two elongated structures, with the measured proper mo- 
tions aligned along the structures' elongation axes. Using 
the same nomenclature of Paper I, we refer to such struc- 
tures as Group I and Group II. Each group consists of two 
(blue- and red-shifted) clusters of features separated by a 
few hundreds of AU for Group I and thousands of AU for 
Group II. 

We detect a compact radio continuum source at 8.4 
and 22 GHz (previously reported at 4.9 and 8.4 GHz) 
located exactly at the center of the Group I maser dis- 
tribution, indicating that the source ionizing the gas is 
also responsible for the excitation of the water masers. 
Comparing our VLA-B observations with previous VLA- 
C and VLA-D data, there is indication that the continuum 
emission is extended on angular scales > 4" and is opti- 
cally thin at frequencies larger than 8.4 GHz. The most 
likely interpretation is in terms of a young HII region. 

The measured transverse velocities of maser features 
indicate that the two clusters of Group II are moving 
away from each other along the cluster-connecting line. 
The spatial and velocity distribution suggests that the wa- 
ter masers of Group II trace a collimated outflow. Since 
the axis of the larger-scale HCO + - SiO outflow is almost 
parallel with the direction of elongation and motion of the 
Group II maser features, these latter likely trace the in- 
ner portion of the outflow and originate in dense clumps 
displaced along the outflow axis. 

Group I maser features individuate a clump of dense 
gas around the expected YSO position (pinpointed by the 
VLA continuum) and the measured proper motions ap- 
pear to trace expansion. If a single YSO were responsi- 
ble for the excitation of both maser Groups, the proposed 
disk- wind scenario might explain the observed flow motion 
of both Group I masers, tracing the dense material emerg- 
ing from the disk atmosphere, and the Group II maser jet, 
corresponding to the portion of the wind collimated along 
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the disk-axis at larger distances from the YSO. Since the 
axes joining the two clusters of maser Groups I and II are 
not intersecting at the position of the continuum location, 
we favour the scenario where two distinct massive YSOs 
are driving the two maser outflows. Presently, the avail- 
able observational data, from radio to NIR wavelengths, 
do not permit to identify unambiguously more than a sin- 
gle massive YSO in the region of the detected maser emis- 
sion. Our detection of the maser outflows can drive future 
studies aiming to localize the powering YSOs. 

Appendix A: Kinematical model formalism 

In order to investigate the kinematical structure traced 
by Group I water masers, we fitted the positions and 
the 3D (line-of-sight + transverse) velocities of maser fea- 
tures with two simple alternative kinematical models: a 
Keplerian disk and a conical outflow. The model fitting 
was performed using both measured absolute and relative 
velocities. 

A.l. Fit of the absolute velocities 

The best-fit solution is found minimising the \ 2 given by 
the sum of the square differences between the model and 
measured, absolute, 3D velocities: 



model predicted and measured relative (to the reference 
maser feature) velocities: 



nspot 
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(A.1) 



where the lowercase v denotes the velocity component 
computed with the model and the uppercase V is the mea- 
sured velocity component, with the corresponding uncer- 
tainty given by AV. In the first term of the sum, only the 
line-of-sight velocities (along the "z" axis) of the identi- 
fied features (indicated with the index i) are considered. 
In order to compare with the model, the measured line-of- 
sight velocities have been first corrected by the peculiar 
motion of the SF R (V^innH = —4.4 km s _1 , as deduced 
by CO emission - ISnell et al.lll988^ . implicitly assuming 
the LSR velocity of the (proto-)star to be the same of the 
surrounding molecular cloud. The uncertainty of the line- 
of-sight velocities is taken equal to the spectral FWHM of 
the maser line, J$ 1 km s . The second term of the sum 
includes the features (indicated with the index k) with 
measured transversal velocities, whose components along 
the R.A. ("x") and DEC ("y") axes are denoted with the 
index j. In the model, the transversal velocity of the YSO 
is taken to be zero, so that we compared directly the model 
and the measured absolute transversal velocities. 



A. 2. Fit of the relative velocities 

Fitting absolute velocities has two disadvantages: 1) the 
absolute motions have large errors (see Col. 8 of Table[3|); 
2) the model assumption that the YSO does not move 
across the plane of the sky may be wrong. In order to 
overcome these limitations, the fit was repeated comparing 
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where the are the measured relative proper mo- 

tions, AIA fe are the associated uncertainties (within a few 
km s _1 in all cases) and the index r indicates the reference 
maser feature. Comparing the model and observed relative 
velocities has the advantage not to require to include as a 
further model free parameter the velocity of the YSO. 
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pointing to us the presence of previous observations indicat- 
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